Cold Atomic Collisions: Coherent Control of Penning and Associative Ionization 
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Coherent Control techniques are computationally applied to cold (1 mK < T < 1 K) and ultracold 
(T < 1 /itK) Ne*(3s, 3 P 2 ) + Ar( 1 S'o) collisions. We show that by using various initial superpositions 
of the Ne*(3s, 3 P2) M = {—2,-1,0,1,2} Zeeman sub-levels it is possible to reduce the Penning 
Ionization (PI) and Associative Ionization (AI) cross sections by as much as four orders of magnitude. 
It is also possible to drastically change the ratio of these two processes. The results are based on 
combining, within the "Rotating Atom Approximation", empirical and ab-initio ionization- widths. 

PACS numbers: 



Cold and ultracold atomic processes present a new lab- 
oratory environment in which to explore and utilize the 
quantum nature of matter. In this letter we show that 
the significance of quantum effects in such systems per- 
mit unprecedented control over collisional processes. In 
particular, we consider the theory and computational im- 
plementation of the coherent control0,0,U0] of absolute 
and relative cross sections in the collision of mctastablc 
atoms A* and ground state target atoms B. Such colli- 
sions can result in two main channels: (1) the ionization 
of the target atom and the de-excitation of the metastable 
species, i.e., Penning Ionization (PI) Q, or (2) Asso- 
ciative Ionization (AI), wherein the colliding partners 
form an ionic dimer while emitting an energetic electron. 
Schematically, 



(AI) AB H 



A* + B -> A + £ H 



(PI) (1) 



As an example, we consider the coherent control of PI 
and AI resulting from collisions between Ne*(3s, 3 P2) and 
Ar( 1 S'o) in the cold and ultracold regimes. Amongst 
other observations, the suppression of these processes 
in favor of elastic scattering may well prove useful for 
the production of Bose-Einstein Condensates of excited 
states atoms. 

A wealth of experimental information and theoretical 
calculations on the uncontrolled Ne* + Ar collisions is 
available and thepossibility of control of this system at 
thermal energies^ now exists. As we report below, the 
control achievable in the sub mK regime is far more im- 
pressive. 

The Initial Superposition State: Coherent Control is 
achieved by preparing the colliding pair in an initial su- 
perposition of internal states, such as, 



|V) =e 



iK-RcM+ik-r 



(2) 



M 



are Ne* Zeeman sublevels, with M = {—2,-1,0,1,2} 
being the projection of the Ne* electronic angular mo- 
mentum on the space-fixed quantization axis, olm are 
preparation coefficients, to be optimized to yield a de- 
sired objective, Rcm is the CM coordinate, Rcm = 
(mN e r Nc + mA r r Ar ) /(thn c + mAr) > and r is the inter- 
nuclear separation vector, r = r Nc — r Ar . The (Body- 
Fixed) momenta are given as, K = k Nc + k Ar , k = 
(mArk Ne — TON C k Ar ) /(mNe + «^Ar) ■ Here r Ar and k Ar 
(r Nc and k Nc ) denote the position and momentum of the 
Ar (Ne*) atom in the laboratory frame. Note that the 
fact that the initial superposition state is comprised of 
degenerate M states, and that the collision partners are 
atoms, ensures that the conditions for coherent control^ 
are satisfied. 

The rates of the PI and AI processes mainly depend 
on A, the body-fixed (BF) projection of the electronic an- 
gular momentum on r, the interatomic axis. It is there- 
fore necessary to express the \4>N e *) states in terms of 
the |<^Ne*) states. We adopt the "Rotating Atom 
Approximation" according to which the axis of quan- 
tization of the electrons faithfully follows the intcrnuclcar 
separation vector. This establishes a 1:1 correspondence 
between the M values and the A values as the atoms 
approach one another. Hence, the (even parity) linear 
combination in the BF frame is written as, 



I^Ar}^ 



\<f>Nc*) a n, 

0=0 



(3) 



where |"0Ar) is the initial state of the Ar atom and \4>^ c »), 



where SI = |A|, and (due to the assumed even parity) 
an = (cm + a- M ). 

Scattering Theory: The basic formulae for our pur- 
poses are found in Refs. 0, 0, 0] giving the scattering 
amplitudes for PI and AI based on O'Malley's theory of 
dissociative attachment 0, 0|. Prior to the collision 
the internuclear momentum vector k has magnitude k 
and direction k. After the collision its magnitude is kf 



2 



and its direction is k/. Asymptotically the Penning elec- 
tron departs along the k e direction with energy e. The 
energy of the emitted electron is related to the collisional 
energy E and the energy of the nuclei after the collision 
E' by the conservation of energy E + Eo = e + E' , with 
co = E„ — IE being the difference between the excitation 
energy of the metastable Ne* atom and the ionization 
energy of the target Ar atom. 



The scattering amplitude, which is exact within the 
Born-Oppcnheimer approximation, is given by IfSLllflj 

/(k/, £ ,k £ ; k) = - 2 J^ (^) 1/2 (^\v e ^m, 

(4) 

where M r is the reduced mass of the nuclei and p e is 
the density of electronic continuous states. V'd( r ) is the 
incoming wave function calculated on the optical poten- 
tial V*{r) — hT{r) and tj} s {r) describes the system on the 
exit channel V+(r). The electronic part is completely in- 
cluded in V £ ^ (r), which is the probability amplitude for 

the emission of an electron with e, and k F . 



Partial wave expansions of ips(r), V e ,i(r), and ipd(r), 
and the evaluation of the integral over r will give, for the 
special case when the space fixed z-axis is along fc, 
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for AI 5j. The V^j(r) matrix elements are related to 
the ionization widths, T n (r) as H H, 13, V$(r) 
a t [T n {r)/2ii] 1 / 2 . The Tsiare obtained as in Ref. 0], 
using methods in [H E EE E3 • 

The differential cross section for PI is obtained by 
squaring the scattering amplitude Eq. (J5J). In the ro- 
tating atom approximation we obtain 



<jg(k/,e,k e ; fc,a n ) = 



^ao/ 9 n(k/,£,k £ ;fc) 
n 



(10) 



where q indicates the exit channel. There are three possi- 
ble exit channels characterized by the electronic state of 
the products: X 2 S^ 2 , A\ 2 Hs/2, and A 2 2 Si/2- The en- 
trance channel optical potential are connected with each 
of the exit channels using Morgner's T splitting [l8lll^ |. 
We use optical potentials for the entrance channel de- 
rived directly from experiment [Til fl5| . Since data is 
available only for scattering experiments at thermal en- 
ergies (above IK), the optical potentials do not include 
effects associated with very slowly moving atoms (e.g., 
hypcrfinc interactions). 

The sum Eq. I|1U|I can be expanded to give 



0"g(k/,£,k e ;fc,an) = y^^a^,a n a q g, n 



(11) 



(5) 



with the partial-wave S-matrix in terms of the phase 
tfj l d and ip l 



shifts S l and Sh of the radial partial wave components 



r, 2 



M+tfH^VetM). (6) 



where V et {r) w a e [T{r) /2n] 1/2 . 



For crossed beams, in the rotating atom approxima- 
tion, we find that the scattering amplitude for a linear 
superposition can be written as 

f a (kf,e,k £ ,k) = y^an/n(fc/,e,k E , fc), (7) 



with fa(kf, e, e; k, f2) given by equation JSJ and S\, t (e) 
replaced by 

Sl;f(e) = -2^^^ sl f\4\V^ l d ), (8) 



with <JqQ,'VL = fqti'fqn- Since we are interested in the total 
ionization cross section we sum over all the exit channels, 
and integrate over the solid angles k e and k/, and over 
the emitted electron energy e. It is easy to see from the 
form of the scattering amplitude that the integration over 
the solid angles yields a Kronecker delta function, thus 
simplifying the expressions for the fgQ>f q n products to: 

/^/*> = £E( 2l+1 K 2 ' /+1 )(o o o) sfi'b^sftb' 

(12) 

The PI cross section, cr PI (a fi ), obtained after integrating 
over the energy of the emitted electron and summing over 
all exit channels, is, 



-Pi 



(an) 



E* pi 



(13) 



As a first example we examine coherent control ob- 
tained using only two f2(= 0, 1) states [similar results 
were obtained for the 0(= 0, 2) pair, as shown in the 
tables below], for which, 



.pi 



(a n ) = \a 



2 PI 



lad 2 erf 1 



2M(a*a 1 a*l) 1 



(14) 
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where a& = oj^. 

Similar expressions for AI are obtained by summing 
over the exit channels and bound states, 

<j m = |a | 2 a AI + \a 1 \ 2 <7 J t l + 2M(a* a 1 a$). (15) 

Note the crucial interference term, dependent on the 
magnitude and phase of the a j , which allows control over 
the cross sections by varying these coefficients. 

Computational Results: Although results are reported 
for collisions at temperatures up to 1 K, our main fo- 
cus is on cold collisions at a temperature of 1 mK and on 
ultracold collisions at 1 /iK. At the low temperatures con- 
sidered, the PI or AI cross sections are very large since 
the two atoms are in the vicinity of one another for an 
extended period of time. Further, in the ultracold case 
only the s partial wave contributes, with three angular 
momentum states contributing in the cold case. 

For these energies the relative velocities between the 
collisional pair are ps 1 m/s and f=a 0.006 m/s respectively. 
These relative velocities are experimentally attainable us- 
ing laser cooling and manipulation techniques. For exam- 
ple, the atoms can be cooled and trapped in a 3D optical 
lattice, then adiabatically accelerated along a single axis 
|2c| . This setup can reach velocities of up to a few meters 
per second, and a kinetic energy spread of 150-200 nK 
around the central beam velocity. The present collisional 
scenario would require two 3D lattice setups in order to 
control both scattering particles. The internal state su- 
perposition can be prepared |2l| after cooling while the 
atoms are trapped in the lattice, for this two techniques 
are possible: stimulated Raman adiabatic passage (ST I- 
RAP) [HI or coherent population trapping (CPT) |2^ . 

Consider first control results across a broad spectrum 
of energies up to temperatures of 1 K, Figure ^ shows 
the cross sections for O = and = 1 scattering as a 
function of energy. Also shown are the maximum and 
minimum controlled cross sections (optimized over the 
a, coefficients) at each energy for the O = plus = 1 
linear combination. Several resonances [2^ are evident, 
since the collision energy is very close to the dissociation 
threshold for the Ne*-Ar quasimolecule. Noteworthy is 
that control is extensive, with enhancement and suppres- 
sion of both cross sections being possible at both resonant 
and nonresonant energies. 

Table H] presents numerical results for the cold collision 
(1 mK) case. We see that it is possible to actively change 
the AI and PI cross sections by as much as four orders of 
magnitude for the O = 0, 1 linear combination and three 
orders of magnitude for the O = 0,2 linear combination. 
For both linear combinations the position of the minima 
and maxima for cr PI and <r AI occur at close points in the 
parameter space (not shown here). A similar observation 
has been noted in the thermal case (above IK), indicat- 
ing that both PI and AI cross sections can be controlled 
simultaneously 0. Sample results for the control of cr PI 



as a function of ai for the cold collision case are shown 
in Fig. 

Ultracold collisions, where only s waves contribute to 
the process, show even more dramatic behavior. As seen 
in Table [HI active changes of up to four orders of mag- 
nitude, using the O = 0,l superposition states, and up 
to three orders of magnitude, using the O = 0,2 super- 
position states, are possible. The AI process can also 
be almost as well controlled. The resulting c AI cross 
sections are shown in Fig. \Bp for ultracold collisions as 
function of the a;. Note that in all cases, the maxima 
and minima in the control plots (Figs. H an dEJ are well 
separated, making the experiment less sensitive to the 
control parameters. 

In summary, we have shown the possibility of a huge 
range of control of the PI and AI cross sections in Ne* + Ar 
cold and ultracold collisions. Control is achieved by ini- 
tiating the collision in a judiciously chosen superposi- 
tion of Ne* quantum states. Such states can be read- 
ily made using new STIRAP techniques j^]. Our results 
show a wide range of controllability for both PI and AI. 
For PI the minimum of the cross section is found to be 
orders of magnitude smaller than the incoherent mixture 
of = 0, 1 or O = 0, 2. In the AI reaction the effects are 
even more dramatic, showing a minimum four orders of 
magnitude smaller than that of the incoherent mixture 
O = 0,l and three orders of magnitude smaller than the 
O = 0,2 mixture. 
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TABLE I: Cross section for cold collision at T = 1 mK. Rows 
labelled "fi = 0, 1" and U Q = 0, 2" show the minimum and 
maximum of the a PI and a AI , obtained by varying the a; for 
the indicated superposition. 



n 


(T PI (A 2 ) 


er AI (A 2 ) 





74.68 


346.91 


i 


64.90 


306.25 


2 


13.75 


87.01 


0,1 


1.27xl0- 2 - 139.57 


3 x 10~ 2 - 653.13 


0,2 


0.63 - 87.80 


0.60 - 433.32 



TABLE II: As in Table I, but for ultra cold collision at T = 
ljuK. 





cr PI (A 2 ) 


<t ai (A 2 ) 





1357.32 


7056.92 


i 


1174.19 


6199.40 


2 


244.27 


1728.65 


0,1 


0.128 - 2531.38 


3.03 - 13256.30 


0,2 


9.88 - 1591.71 


2.28 - 8783.30 



250 





FIG. 1: PI and AI cross sections for SI = 0, 1 linear combina- 
tion, fl — -solid; SI = 1 -dashed; maximum and minimum 
for the linear combination of f2 = 0, 1 -dotted. 
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FIG. 2: Coherent Control contours for (a) PI for SI = 0, 1 in 
cold collisions at T = 1 mK; (b) AI for SI = 0, 2 in ultra cold 
collisions at T = 1 /xK. The parameters rj and A£ are defined 
via: an = sinrje^ n and a n / = cos 77 e 1 ^"' , with f2, Q' — 0, 1, 2. 



